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ABSTRACT 



Thermospheric airglow emissions in the wavelength range of 1800 A to 3400A and 
altitude range of 100 to 330 km were identified. These emissions compose the preliminary 
data analysis for the Naval Post Graduate School MUSTANG rocket experiment scheduled 
for launch from White Sands, New Mexico in February of 1990. Identification of 
emissions in the wavelength range of the MUSTANG experiment was accomplished using 
experimental results to support calculated emission intensities. All emission bands that will 
contribute to the experimental spectra, as well as those features that will be distinguishable, 
have been tentatively identified. These results confirm the contention that the instrument 
wavelength range allows observation of the major midultraviolet airglow emission bands. 
The extended wavelength and altitude range for this experiment allows a unique 
experimental opportunity to observe all major neutral atmospheric constituents in one 
spectrum. The primary emission band systems that will be observed are the N2 Lyman- 
Birge -Hopfield bands, the O2 Schumann-Runge bands, the NO Gamma bands and the O 
2972A line. Secondary contributions will be observed from the N2 Vegard-Kaplan and 
2 nd Positive bands, the O2 Herzberg bands, and the NO Epsilon and Delta bands. These 
combined spectra will insure that the MUSTANG experimental goal of measuring N2 and 
NO densities and determining EUV and photoelectron fluxes can be met. 
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I. INTRODUCTION 



There are many factors motivating the investigation of ionospheric photochemistry. 
The Department of Defense has several areas of interest in this field. A 1986 Joint Chiefs 
of Staff memorandum identified fifty critical geophysical parameters that need to be 
determined. The determination of electron density was ranked fifth on the list (MJCS, 
1986). The DoD's concern in furthering the understanding of the ionosphere is based on 
three known needs: 

- Measure the ionosphere for over the horizon radar applications. 

- Measure the ionosphere for uni-directional communications. 

- SDI beam weapons require identification of low absorption windows for 

propagation in the upper atmosphere. 

Additionally, there is a purely scientific interest in fully identifying physical processes in 
the upper atmosphere. 

These motivating factors have led to the Naval Post Graduate School rocket 
experiment called MUSTANG. The mission of this experiment is to investigate 
ionospheric photochemistry. Three specific experimental goals have been established. 
These are: 

- measure NO density in the E-region ionosphere. 

- measure N 2 density in the E-, FI, F2 regions. 

- infer solar EUV flux and photoelectron flux. 
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This is a joint rocket experiment containing two ultra-violet spectographs. The 
Naval Post Graduate School is responsible for a spectograph measuring ultraviolet 
radiation from 1800 angstroms (A) down to 3400A. The Naval Research Laboratory is 
responsible for a spectograph measuring from 500A down to 1500A. The launch vehicle is 
a Black-Brandt sounding rocket, which is scheduled for a February, 1990 launch date at 
White Sands, New Mexico. NASA is responsible for the launch. Observations will be 
made from 100 to 330 km in altitude. 

A. THESIS GOALS 

The primary objective for this thesis is to identify emission band systems with 
features between 1800 and 3400 A in wavelength, that are observable from 100 to 330 km 
in altitude. The accomplishment of this thesis objective will assist in verifying that the 
MUSTANG experimental goals are attainable. 

B. THESIS OUTLINE 

Chapters II through V lead to the identification of those atmospheric constituents 
with emission features in the 1800-3400 A range. Chapter II describes the experimental 
environment which is the ionosphere between 100 and 400 km. Chapter III describes the 
general physical process that are the principal methods of transferring solar flux energy to 
create the ionosphere and airglow. Chapter IV describes two facets of molecular 
spectroscopy. The first of these is molecular transitions. The second is synthetic spectra. 
Chapter V identifies which emission band systems for the ionospheric constituents (N 2 , 

O 2 , NO, O) are observable in the experiment. Finally, Chapter VI provides conclusions 
and suggestions for further research. 
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II. THE IONOSPHERE 



The purpose of this chapter is to describe the region in space that is the focus of this 
experiment: the ionosphere. It is helpful to consider the entire atmosphere, before 
proceeding to focus on the ionosphere. The general atmospheric regimes are shown in 
Figure 2.01. Several of these regimes will be discussed in this paper. The first two 
regimes of concern are the homo- and heterospheres. The division of the atmosphere into 
these two sub-spheres is a result of significantly different processes that influence the 
layered structure of the ionosphere. Therefore, the homosphere and heterosphere will be 
discussed initially, followed by a description of the ionosphere, its sources, and its 
sublayers. 

A. THE HOMOSPHERE AND HETEROSPHERE 

The homosphere is the region of the atmosphere where turbulent mixing or eddy 
diffusion is the dominant dynamic mixing process. Atmospheric constituents N 2 , O 2 , Ar, 
He, and H, are uniformly distributed in this region. The homospheric air density falls off 
exponentially with increasing altitude as shown below: 



p (z) = po (2.1) 

The altitude is given by z and H is the atmospheric density scale height. This is defined 
below as 



H BjlL 1<>I 

n m g + T dz 



( 2 . 2 ) 
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where kg is the Boltzman constant, T is atmospheric temperature, g is acceleration due to 
gravity, and m is the mean molecular weight of the medium. This scale height is common 
to all constituents in the homosphere. The expression defining the scale height can be 
simplified by eliminating the second term when the temperature profile becomes isothermal 
and the temperature gradient goes to zero. 




Figure 2.01 Atmospheric Regimes (Tascione, 1988) 
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The turbulent mixing processes associated with the air-earth interface diminish with 
increasing altitude. At around 100 km the dominant mixing mechanisms become molecular 
diffusion rather than eddy diffusion. As shown in Figure 2.01 the boundary is known as 
the turbopause. Atmospheric constituents begin to separate out based on molecular 
density, much as oil and water separate. The number density for each constituent falls off 
exponentially in proportion to its scale height. The density scale height is no longer based 
on mean molecular weight, but rather each individual molecular weight. This region is 
called the heterosphere. 

B. THE IONOSPHERIC ENVIRONMENT 

With increasing altitude, molecular diffusion rather than eddy diffusion becomes the 
dominant mixing process. This process combines with an increasing solar flux, and a 
decreasing atmospheric density to produce ions that have a sufficiently long lifetime to 
support the existence of a continuous plasma. This plasma is known as the ionosphere. 

The general relationship between incident solar flux and atmospheric gas concentration is 
shown in Figure 2.02. As is shown in this figure the lower limit of the ionosphere is 
typically 50 to 70 km (Tascione, 1988). The relatively wide range is primarily due to 
fluctuations in solar flux. The bold curve in Figure 2.02 represents the number density of 
ions over altitude The ionizing radiation is primarily the extreme ultra-violet and x-ray 
portion of the solar spectrum. The general shape of the ionospheric region indicates that 
ion density is negligible at very high altitudes. It then begins increasing in number density 
to a peak and decreases to a negligible intensity by an altitude of 50 km. 
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Figure 2.02 Generalized Ionosphere 
The fluctuation in the upper and lower limits of the ionosphere, as well as 
fluctuations in the ion density vary in proportion to the level of solar activity. The vast 
majority of solar flux may be treated as a black body spectrum, based on the solar surface 
temperature, of about 5700 degrees kelvin. However, the high energy end of the 
spectrum, in the ultraviolet and x-ray regions, as well as the low energy end at radio 
frequencies exhibits significant fluctuations (Torr et al, 1979 for solar flux data over the 
solar cycle). The composition and vertical limits of the ionospheric layers are significantly 
influenced by this varying solar flux. The variable portion of the solar flux varies in 
predictable cycles. This solar cycle has a period of approximately eleven years. An 
accepted indication of solar activity is the measurement of a solar radio emission whose 
wavelength is 10.7cm. This measurement is designated the F10.7 and is taken at Ottawa, 
Canada. Both the F10.7 number and the sunspot number are shown in Figure 2.03. As 
the figure indicates, the F10.7 number tracks closely with the sunspot number. 
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Figure 2.03 F10.7 and Sunspot Index(NRL, 1989) 

The variable solar flux results in atmospheric temperature fluctuations, especially in 
the upper region of the atmosphere known as the thermosphere. The thermosphere is 
primarily heated by atomic oxygen absorption of extreme ultra-violet radiation (Tascione, 
1988). Because, the extreme ultra-violet flux varies widely, the temperature of the 
thermosphere also varies. As is shown in Figure 2.04 this results in thermospheric 
maximum temperature ranges from 500 to almost 2000 degrees Kelvin. The predicted level 
of solar activity allows an atmospheric temperature profile to be estimated for a given 
location and date. An atmospheric temperature profile has been prepared for the 
MUSTANG experiment in Figure 2.05. The profile was obtained using the neutral 
atmospheric model MSIS-83 (Hedin, 1983). The model assumes a launch site latitude 33° 
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N, and longitude 105° W at White Sands, New Mexico. It also uses the predicted FI 0.7 
value for an early February 1990 launch date of 215. The model also considers the launch 
time of day, which will be 1000 hours. The solar zenith angle is 72.47 degrees. Active 
solar conditions are used in the model. 




KINETIC TEMPERATURE l°K) 



Figure 2.04 Atmospheric Temperature Ranges (Tascione 1988) 
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Figure 2.05 Projected Launch Date Temperature Profile (Hedin, 1983) 

Solar radiation of sufficient energy to ionize atmospheric constituents is the primary 
ionizing mechanism throughout the ionosphere. The far and extreme ultraviolet, portion of 
the solar spectrum contain enough energy to ionize the principal constituents of the 
ionosphere. The minimum ionizing energy requirements for the constituents are 
determined by the ionizing work function of each species. These work functions are given 
in Table 2.1 for the primary atmospheric inhabitants 
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TABLE 2.01 



IONIZING WO RK FUNCTIONS FOR NE UTRAL PARTICLES 

ITascione, 1988) 

km KE teY) UA1 



n 2 


15.58 


795 


N 


14.50 


851 


O 


13.20 


908 


O 2 


12.07 


1027 


NO 


9.26 


1339 



Table 2.02 indicates that the ionizing radiation requirement is at the highly variable, x-ray 
and extreme ultra-violet portion of the solar spectrum. The principal ionizing portion of the 
solar spectrum is shown in Figure 2.06. It is interesting to note that the strong Lyman 
Alpha solar emission at 1215A will only meet the ionization energy requirement for the 
brelatively scarce NO molecule which results in its increased importance in especially the D 
region of the ionosphere. Photoionization is a direct energy absorption. The 
photoionization reaction is shown below for molecular oxygen: 

02 + h\) => 02 + + e 
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Figure 2.06 Solar Radiation 250-1250 A (Banks & Kockarts, 1973) 



This is the initial ionization process for each of the major atmospheric constituents. The 
kinetic energy of the electrons produced may be quite significant. These electrons are 
known as, photoelectrons. They will participate in one or even several additional 
collisions, before finally reaching thermal energies of around .025 electron volts (Banks 
and Kockarts, 1973). 

C. THE IONOSPHERIC LAYERS 

A closer examination of the ionospheric zone shown in Figure 2.02 reveals that it is 
not a smooth curve. The vertical structure of the ionosphere is changing constantly. It also 
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varies with time of day, solar activity, as well as with latitude. Despite fluctuating 
dramatically over time, and location, the ionosphere is always present. Furthermore, the 
ionospheric plasma has four distinct regions which are primarily the result of three 
competing factors: 

- the atmospheric absorption of solar flux varies with altitude 

- recombination of ions is density dependent 

- atmospheric composition varies with altitude. 

The ionospheric layers are designated the D, E, FI, and F2 region. The layer 
boundaries are not distinct. Although they are all present during the day, some do not 
survive throughout the night. Their general peak intensities correspond to the altitudes 
shown in Table 2.02. 

TABLE 2,02 

IONOSPHERIC REGIO NS ; PEAK HEIGHTS ( Km) 

D -90 
E -110 
FI - 200 
F2 - 300 

These layers, are called Chapman layers after Sidney Chapman, who first modeled them in 
1931 (Chamberlain, 1987). The general shapes and locations of the Chapman layers are 
shown in Figure 2.07. The solid lines in the figure indicate electron concentrations under 
active solar conditions. The dashed lines indicate solar minimum conditions. 
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DAY/NIGHTTIME ELECTRON CONCENTRATIONS 




Figure 2.07 Chapman Layers (Chamberlain, 1987) 



Before describing each ionospheric region it is necessary to describe the distribution 
of neutral species in the heterosphere, which overlays the E and both F layers of the 
ionosphere. The relative abundance of the neutral particles is distributed differendy in each 
layer. The neutral atmospheric particle distribution is shown in Figure 2.08. Below about 
175 km N 2 is the dominant species present. Above this altitude atomic Oxygen becomes the 
dominant species. O 2 and NO remain relatively minor constituents although as will be 
shown in Chapter V, both have significant roles to play in ionospheric photochemical 
processes. 
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Figure 2.08 Primary Ionospheric Constituents 
(Hedin, 1983 and Clear>. 1986) 

The dominant processes that produce each ionospheric region is described below. 
A more detailed treatment of each region may be found in The Theory of Planetary 
Atmospheres , by Chamberlain (1987). 

1 . The D Region 

The D region extends from 50 to about 90 km, with its peak near 90 km. 
The relatively high number density in this region results in a short mean free path between 
neutral particles, and electrons. The daytime electron density is about 1 .4x 10 4 electrons 
per cubic centimeter (el cm* 3 ). This produces a high recombination coefficient that causes 
the region to disappear at night. Solar Lyman-alpha radiation deposited in the region is the 
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principal ionization source for the NO molecule. Although photoionization is the dominant 
process, other interactions contribute to ion density. These include primarily three-body 
attachment as shown below: 

O 2 + e + M => O 2 ' + M. 

A reaction of secondary significance is radiative attachment such as the example shown 
below: 

O 2 + e => 02 ‘ + hi). 

Principal loss mechanisms include photodetachment, which is the inverse of radiative 
attachment; collisional detachment, the inverse of three-body attachment; and mutual 
neutralization, such as X + + Y* =* X + Y. This last process indicates that the D region is 
characterized by a significant number of negative ions as well as electrons. 

2 . The E region 

The E region extends from 90 to about 130 km. The nominal peak height is 
1 10 km. The daytime density of 1.5 xlO 5 (el cm -3 ) decays down to a night time density of 
lxlO 4 (el cm* 3 ) or less. The decreasing recombination coefficient allows its nighttime 
survival. Ionization primarily is the result of photoionization by solar extreme ultraviolet 
and soft x-rays. The ionization of atomic oxygen becomes significant. Negative ions are 
no longer significant contributors in the region. The principal ion loss mechanism is 
dissociative recombination such as: 

NO+ + e => N +0 

3. FI region 

The FI region extends from 140 to 200 km. Its daytime electron density is 
greater than the E region (2.5x10^ el cm' 3 vs. 1.5x10^ el cm* 3 ) However, the principal 
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ions formed are 0 + and NO + which require an ionizing solar flux from the Lyman 
continuum as well as He emission lines. These sources disappear with sunset. Due to the 
relatively short ion lifetimes, the FI layer disappears quickly at night. 0 + readily charge 
exchanges with O 2 and NO. These exchanges are given below: 

0+ + O 2 =s- O + 02 + 

0+ + N 2 => N0+ +N. 

Once in these molecular forms the principal loss method is once again dissociative 
recombination. 

4 . F2 region 

Beginning at 160 km, the F2 region extends all the way up and beyond 600 
km, with a nominal peak density at about 300 km. Its day time density is around 10^ (el 
cm‘3). The night time density drops by about one order of magnitude, to about 10^ (el cm" 
3 ). In the F2 region three significant factors adjust the ionization profile. First, the 
atmosphere becomes optically thin to most ionizing radiation; which adjusts the ion 
production equation. Second, ion-atom interchange becomes the dominant loss process 
instead of dissociative recombination. This results in the loss term diminishing as 
atmospheric density decreases. Ion production is from the same processes as the FI layer. 
Atomic oxygen is the principal constituent and ion source in the region. However, since 
the loss term is diminishing with altitude the ion number density would increase indefinitely 
with height except for the effect of the third factor, ambipolar diffusion. The F2 region 
reaches the point where the chemical lifetime becomes longer than the diffusion lifetime. 

As a result positive ions diffuse downwards out of the region. This limits ion density 
towards the upper end of the region. 
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III. THE SOURCES OF THE ATMOSPHERIC AIRGLOW 



A. BACKGROUND 

The MUSTANG experimental objective to investigate ionospheric photochemistry 
is best accomplished by observing radiative emissions in the earth's airglow. Airglow 
emissions are omni-directional. An instrument's viewing angle is oriented to observe 
emission features without being overwhelmed by the solar flux. Specific emission features 
can be observed to identify the atomic or molecular transitions and therefore the molecular 
density of the species. Airglow emissions can then be converted to density profiles. These 
number densities can then be used to infer the electron and ion densities. 

The sum of all radiative emissions from the extreme ultraviolet to the infrared 
comprises the atmospheric airglow. The airglow is optical radiation that arises from several 
types of interactions. Airglow emissions have three general sources: 

- direct sunlight scattering 

- electron impact excitation 

- photochemistry of neutral constituents. 

The primary source is direct sunlight scattering, or photoexcitation. Electron impact 
excitation is of secondary importance for allowed transitions, but it is the primary excitation 
mechanism for forbidden transitions. Some radiation is the result of a combination of these 
factors. 

Solar ultraviolet radiation is primarily responsible for the airglow to be observed in 
the MUSTANG experiment. Figure 3.01 shows a typical nadir spectrum of the earth's 
airglow from 500 A to 4000 A. The primary EUV and FUV emission features are shown in 
this figure.The strong O 2 and O 3 absorptive bands shown in the figure result in the 
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suppressed airglow at wavelengths corresponding to these absorption bands. The figure 
also illustrates the emission bands to be observed by the Naval Research Laboratory 
instrument . The strong Lyman Alpha, and atomic Oxygen emission features are clearly 
visible between 500 and 1500A. Conversely the NO and N2 features in in the 2000 to 
3000 A range are not visible in a nadir view. This is an illustration of why the experiment 
must look across the limb to observe emissions in the 1800 to 3400 A region. 
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Figure 3.01 Earth Ultraviolet Airglow and Absorption/Emission Regions 
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The airglow is distinguished from another emission phenomenon, the aurora, by 
the energy source of the aurora. The energy sources for the aurora are solar electrons and 
protons, as opposed to solar flux as the primary source of the airglow. This results in the 
auroras confinement to polar regions and its sporadic occurrence. Conversely the airglow 
occurs continuously, over the entire sunlit portion of the atmosphere. 

A. DIRECT SUNLIGHT SCATTERING 

The primary process that results in the airglow is direct sunlight scattering, or 
photoexcitation. This process is similar to photoionization, the major ionization source. 
Photoexcitation occurs when the incoming photon energy matches the quantum energy 
levels required for a given atom or molecule. The photon will then be absorbed and leave 
the particle in an excited state. The normal lifetime of a particle in an excited state is very 
short; about 10' 8 seconds (Eisberg, 1985). 

The particle in the excited state can return to its desired ground state by emitting one 
or several photons. Photoexcitation that leads to the emission of photons is known as 
fluorescence. The entire process is known as fluorescent scattering. A special case of 
fluorescent scattering is resonant scattering. In the case of resonant scattering, the 
excitation and emission paths are the inverse of each other. A single photon, of the same 
wavelength as the excitation photon, is reradiated. The particle return path in general, may 
be along a direct or an indirect path. Multiple photons are emitted, for the general case 
where intermediate states lie between the excited state and ground. Only electronic 
excitation levels are available to atoms. Molecular excitation may result in electronic, 
vibrational, or rotational transitions, or combinations of these three. Energy requirements 
vary for each type of transition. Electronic transitions require between one and ten electron 
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volts, while vibrational transitions require between one tenth and one electron volt, and 
rotational transitions are a factor of ten less than vibrational transitions. 

Transitions are considered allowed or forbidden depending on their compliance 
with quantum selection rules. The actual selection rules that apply to molecular transitions 
will be discussed further in Chapter IV. A transition that meets all selection rule 
requirements is an allowed transition. A transition that violates a selection rule at one level 
may be forbidden; but may or may not be strictly forbidden. For example a transition 
forbidden as electric dipole radiation may be allowed with a lower probability as a magnetic 
dipole, or electric quadrupole transition. 

An excited state whose only path to the ground state is through a forbidden 
transition is called a metastable state. Usually, the particle can still reach the lower state, 
but the probability of doing so becomes much smaller than if the lower state was an 
allowed state. Lifetimes in the excited state can become relatively long. They can be on the 
order of seconds as is the case for one of the emission bands to be examined in Chapter V. 
Metastable states are much more likely to transfer energy as a result of a collision with 
another particle prior to emission. This process is called quenching. 

It is also of note that molecules in an excited state may dissociate, or predissociate, 
rather than emitting a photon. An absorption leading to an excited state above dissociation 
or predissociation energy, will dissociate rather than emit a photon. These processes are 
best shown graphically. Figure 3.02 shows two simplified potential energy diagrams for 
the O 2 and NO molecules. As shown in the figure, an O 2 molecule excited to the A 3 E + U 
state at a vibrational level above the dissociation threshold energy will dissociate to the 
0( 3 P) + 0(3p) atoms. This may be viewed as a path to climb out of the electronic potential 
well. Conversely, a molecule of NO excited to a vibrational level above the 
predissociation energy of the C^fl state may predissociate along the a 4 Il path to become 
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N( 4 S°) + O(-’P) atoms. This may be viewed as a hole in the potential energy well. The 
result of both of these processes is that absotptions will not result in emissions. 



0 2 Potential Energy Diagram 




NO Potential Energy Diagram 




Figure 3.02 Dissociation and Predissociation Paths (Chamberlain, 1987) 
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B. ELECTRON IMPACT EXCITATION 



The primary excitation source for atoms or molecules is solar radiation. However, 
photoexcitation is forbidden for some of the constituents in the upper atmosphere. In this 
case, electron impact can cause the excitation of the particle to an excited state where it can 
than emit a photon upon undergoing a forbidden transition. Electron impact scattering is an 
inelastic collisional process. The probability of an emission resulting from an electron 
collision with a neutral particle depends on several factors. Electron impact excitation is 
dependent on the molecular number density, electron impact cross-section, and the electron 
flux. Since the electron flux is not monoenergetic, and since the molecular electron 
absorption cross-section is dependent on electron energy, an actual calculation requires an 
integration over available energies. This will be done for several forbidden emissions in 
Chapter V. 

C. PHOTOCHEMICAL REACTIONS 

A third relatively minor source of airglow emissions is radiation resulting from the 
photochemistry of neutral constituents. Photochemical processes are initiated in the same 
manner as photoabsorptive processes. A photon is absorbed leaving the absorbing 
molecule in an excited state. This excited state could be an excited rotational or vibrational 
level or an excited electronic state. It could also be a dissociated or ionized state. In this 
excited state the molecule may enter into secondary chemical reactions that would not be 
available from the ground state. Upon completion of the chemical interaction the new 
species emits a photon. This process is called chemiluminescence. 

Photochemical processes are complicated because they are non-local 
thermodynamic equilibrium processes. Each rate of reaction must be calculated separately. 
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Addressing each photochemical process is beyond the scope of this thesis. 
Chemiluminescence is an insignificant contributor to the day airglow. It becomes important 
at night when other excitation sources are not present. 



23 



IV MOLECULAR SECTROSCOPY 



Chapter IV addresses molecular transitions, and synthetic spectra . These are two 
broad regions that comprise molecular spectroscopy. The vast majority of emission 
processes to be observed in the MUSTANG experiment are a result of molecular 
transitions. Furthermore, this chapter is exclusively concerned with diatomic molecules. 
Diatomic molecules are the dominant type found in the atmosphere. Molecules other than 
diatomic introduce coupling schemes beyond the scope of this thesis. More specifically, 
emissions will primarily occur as a result of vibrational transitions between molecular 
electronic states. This type of emission process is one of several possible types of diatomic 
molecular transitions. Molecular transitions occur between rotational and vibrational levels 
as well as between electronic molecular states. It may also be between some combination 
of these three. Each type of transition is characterized by quantum energy jumps rather 
than a continuous emission spectrum. All three types of molecular transitions will be 
addressed, with emphasis on vibrational transitions between electronic states. Additionally 
the compilation of synthetic spectra will be addressed. A more detailed account of 
molecular transitions is provided by Herzberg, in his definitive volume, Spectra of 
Diatomic Molecules (Herzberg, 1950). 

A. MOLECULAR TRANSITIONS 

This section is concerned with rotational and vibrational levels, as well as electronic 
states. It is also concerned with transitions between these states. These transitions range in 
energy from the least energetic rotational transition, to the most energetic electronic 
transition. They are in effect a nested sequence of transitions. A pure rotational transition 
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can occur between two rotational levels in the same vibrational level and electronic state. 
However, as vibrational and electronic transitions come into the picture, subordinate level 
transitions must be taken into account. 

1. Rotational Levels and Transitions 

A diatomic molecule as shown in Figure 4.01, can be represented as a 
dumbbell. The two masses mi and m 2 are connected by a rigid rod of length r with the 
center of mass at s. The two masses represent the atomic nuclei, where the vast bulk of the 
atomic mass is concentrated. As opposed to an atom, the molecule can undergo two forms 
of motion. It can rotate around the center of mass axis and the two nuclei can vibrate 
relative to one another along the intemuclear axis. The two masses are not in fact rigidly 
bound. This perturbation to the rigid rotator and vibrator will be addressed later in this 
chapter. At this point it is adequate to note that the intemuclear axis will vary with the 
speed of rotation, or strength of oscillation. 




Dumbbell Model of a Diatomic Molecule 



Figure 4.01 Dumbbell Model (Herzberg, 1950) 



The simplest form of motion the two components can undergo is rotation 
around the center of mass axis S. Quantum theory requires the velocity of this rotation to 
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take on discrete values. Isolating the molecular motion to pure rotation results in the 
system's angular momentum as 



P = Ico . (4.1) 



The classical kinetic energy of rotation would be E = ^ I 0)2. The classical energy becomes 
the expression 



E 



P2 

” 21 * 



(4.2) 



P is from equation 4. 1 , and the moment of inertia I, is dependent on the reduced mass. 

The solution to the Schrodinger equation yields available energies of the form shown 

E = illiLL±Jl . (4.3) 

8 7t 2 I 

In Equation 4.3, h is Planck’s constant, I remains the moment of inertia, and the rotational 
quantum number J is introduced. J is limited to positive integer values. The energy of the 
system E, will increase quadratically with increasing J as is indicated in Equation 4.3. The 
angular momemtum of the system in some quantum state J becomes 

P= — V J(J + 1) * — J . (4.4) 

2k 2k 

Classical electrodynamics requires a changing dipole moment to produce an emission. Any 
diatomic molecule composed of two heteronuclear atoms will possess a permanent dipole. 
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Diatomic molecules with a permanent dipole moment will produce an emission when 
rotated. The quantum mechanical system indicates that the system will have a series of 
discrete energy levels whose energy increases quadratically as J increases. A transition 
between two of these available energy states will result in a discrete rotational emission 
feature. The wave number v, is introduced at this time. Wave number is equal to the 
reciprocal of the wavelength (1 A) or frequency divided by the speed of light. Wave 
number carries units of cm* 1 (number of waves per centimeter). The wave number for the 
rigid rotator is given by 



E' 

v "he " he 



( 4 . 5 ) 



The upper and lower energy states are indicated by E' and E", respectively, while h is 
Plancks constant and c is the speed of light. The rotational term F(J), is also introduced, 
and depends on the physical characteristics of the molecule and dependent on J. The 
relationship is 



F(J) = hF = BJ(J + 1} * (4 6) 

The constant B is known as the rotational constant of each molecule and is given by 

B = — ~ . 

8n 2 cl 

Only certain rotational transitions are allowed for a given excited state T. 
The available states are governed by what are known as selection rules. Selection rules for 
the rotational quantum number J, are AJ = ± 1. This requires that J vary by one higher or 
one lower quantum level in either absorption or emission for a rotational transition. 
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2 . 



Vibrational Levels and transitions 



A perturbed, diatomic molecule can also undergo vibratory motion. The 
configuration for the molecule shown in Figure 4.01 is still appropriate, however a second 
mode of motion is now considered. A classical picture would consider the two molecules 
as connected by a spring. Perturbing the system would result in simple harmonic motion. 
The equation of motion is a straight forward result of Hooke's Law: F = -kx, where k is 
the force constant . The vibrational frequency of the molecule is defined as v osc = 

where |X is the reduced mass. For a given vibrational state, the vibrational energy level, 
and the resultant vibrational term G(v) is given by 

G ( v) = ^ = to(u + ±). (4.7) 

The factor co, equal to v 0SC /c, is introduced in equation 4.7, and is a constant for each 

molecule. The vibrational quantum number v is also introduced. The vibrational quantum 

number is restricted to positive integers. It is of note that when v is equal to zero; the 
lowest vibrational state, the vibrational energy equals ^ co, rather than zero. Additionally 

the linear dependence on v, results in equidistant energy levels, for increasing vibrational 
levels. 

Vibrational selection rules, restrict available transitions between higher and 
lower states to unity; Av = ± 1. Homonuclear molecules such as N 2 , have a zero dipole 
moment. Vibrational transitions can only occur as quadrupole transitions for these 
molecules. Quadrupole radiation has a very small probability with respect to dipole 
radiation, and is negligible under normal circumstances. The selection rule Av = ± 1 is 
not rigidly adhered to for asymmetric molecules such as NO. Dipole allowed transitions 
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have a weak probability of emission when Av > ± 1. This is as opposed to the rotational 
case where the selection rule is rigidly adhered to. 

The molecule can undergo simultaneous vibrational, and rotational motion. 
The total energy in the rotating oscillator case is approximately T= G(v) + F(J). The period 
of oscillation is on the order of 1(H* seconds, while the period of rotation is about 10*^ 
seconds. This means that the rotational position changes slowly compared to the oscillating 
motion. For a given vibrational energy level, there are a series of rotational levels 
available. The selection rule remains AJ = ± 1 for a specific vibrational level. The series 
of vibration-rotation transitions is called a band. The AJ = 1 case forms one branch, 
known as the R-branch and the AJ = -1 case forms the P-branch. An illustrative example 
of a vibrational-rotational band system, with P and R branches is shown in Figure 4.02. 



29 



J' 

•10 




c , : 




c t : 




« 1 **,**, “s « 


I | ! 


| | j | j" 


1 




1 1 1 1 ! : i i i 


1 


i i i 



I I I I . I I 




Figure 4.02 Vibrational-Rotational Energy Level Diagram 

(Herzberg, 1950) 

3. Electronic States, Nomenclature, and Transitions 
a. Electronic States 

The nuclei of a molecule repel each other. The molecule is held 
together by what is in effect an entire molecular group of electrons. The electron cloud 
interaction with the nuclei results in a strong axial electric field. The orbital configuration 
of this electron cloud determines the electronic energy of the system. The electrons are 
moving much faster than the nuclei, so the electronic energy corresponds to each temporary 
position of the nuclei. Work must be done against the Coulomb repulsion of the nuclei, 
and work is required to change the electronic energy, in order for the position of the nuclear 
elements to change. There must be some minimum value in order for the electronic state to 
be stable. The shape of the potential energy curves for electronic states indicates whether 
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the state is stable or unstable. Without a minimum in the potential energy curve, the two 
atoms will repel each other at any intemuclear distance. This is an unstable state. Figure 
4.03 illustrates several potential curves for different electron states for a L '12 molecule. As 
is shown in the figure, the , 1 £ + u and 1 Z + g states all possess minimum points, and are 
stable states. The 3£ + u state has no minimum and is an unstable state. 




Figure 4.03 Stable and Unstable Potential Energy Curves 

(Herzberg, 1950) 
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b. Electronic Transitions 

In a similar manner in which purely rotational transitions were 
extended to include vibrational-rotational transitions, vibrational-rotational transitions can 
now be extended to include transitions between electronic states. The total energy of the 
molecule in a given energy state is defined as the sum of the rotational, vibrational and 
electronic components. This definition is given by 

E = E e + E v + E r . (4.8) 

The change in energy from one specific state to another is the sum of the energy changes of 
each component as 



AE = AE e + (E v ' - E v ") + (E r ' - E r "). (4.9) 

For any given change in electronic state, there is a wide range of possible vibrational and 
rotational transitions. Thus, the transitions become a band rather than a single line. Figure 
4.04 illustrates the concept of nested rotational, vibrational and electronic levels. This 
figure illustrates the nesting concept between electronic states. For a given transition 
between two electronic states, there are numerous vibrational transitions possible. Each 
vibrational transition is in fact a band of all possible rotational lines for each particular 
upper and lower vibrational level. The sum of these bands form the band system 
representing emission or absorption features between two electronic states. 
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F, gure 4.04 Rotational, Vibrational, and Electronic Levels 

(Eisberg, 1985) 
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c. Electronic State Nomenclature 

The internal molecular electric field is symmetric about the 
intemuclear axis. The molecular electron's motion in this field results in a total orbital 
angular momentum, L. As shown in Figure 4.05, L is not a constant of the motion as it 
precesses around the axis of symmetry. However, the component of L along the 
intemuclear axis is a constant of the nutation. This component designated A i is a better 
indication of the molecular electronic state than L. The electronic configurations, and 
indeed the labeling of successive electronic states, is analogous to the atomic case. The 
symbol A = I MJ = 0,1,2,... L is used to designate a particular electronic state. 
Furthermore, each electronic state is given a Greek letter to distinguish the molecular 
electronic state from its atomic equivalent. The lower four molecular states are: 

0 - I 

1 - n 

2 - A 

3 - <D. 

It should be noted that for Ml > 0, electronic states will exhibit double degeneracy since for 
a given A, Ml can be positive, or negative. 
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Orbital Angular Momentum Precession 
about the internuclear axis r 

Figure 4.05 Precession about the Internuclear axis (Herzberg, 1950) 

Electron spin in molecules results in a multiplet structure for either a 
particular band or an individual line. Just as in the atomic case, an odd or even number of 
total electrons results in a half integer or integer total spin, respectively. The resultant spin 
S is unaffected by the electrostatic field. It will precess in the presence of a magnetic field. 
When A is greater than zero, the orbital motion of the electrons results in the generation of 
an internal magnetic field in the direction of the internuclear axis. The spin interaction, S, 
with this magnetic field will precess in much the same manner as L. Also in much the same 
manner as L, a constant component Ms will exist for a given spin precession. Ms has 
specific values designated Z, where 

Z = S, S-l, S-2,. . . -S. 

This is a total of 2S+1 possible values. This also indicates the multiplicity of the state. In a 
manner similar to atomic classification, the spin is written as a left superscript. Where A * 
0, the sum A + Z is written as a right subscript. 

Some additional points of classification for electronic states are 
noted at this time. The ground state of a diatomic molecule is labeled with X, while 
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successive excited states, of the same multiplicity, are labeled A, B, C, and so on. 
Similarly, successive excited states where multiplicity differs by two are designated by a, 
b, c, . . . 

The nomenclature being discussed above, is repeatedly used in the 
potential energy diagrams in Chapter V where various emission band systems are 
investigated. Emission band systems are generally named for their discoverer. For 
example the N 2 Vegard-Kaplan bands were discovered by Vegard and Kaplan (JPCRD, 
1977). These bands are designated as a3x u + - X^Zg + . This indicates a transition from 

the first excited electronic state A, to the ground state X for the case where both upper and 
lower states are in the X electronic configuration. Since A = 0 for the X state, the right 
subscript of A + X provides no additional information. In this case the right subscript 
becomes u and g which indicates an odd to even electronic state.The spin goes from X 
equal to three to X equal to one. This violates one of the selection rules that govern 
electronic transitions. These selection rules are discussed in the next section, 
d . Selection Rules 

Selection rules governing electronic transitions are not as simple as 
those for rotational or vibrational levels. The selection rules that govern allowed transitions 
are shown in Table 4.01. Allowed transitions will comply with these selection rules. This 
will result in electric dipole radiation. The AJ selection rule gives rise to three branches 
although more branches may result from a particular transition. Depending on the 
multiplicity of the upper and lower states the three branches are designated R, P and Q. 

The + <=> - selection rule in Table 4.01 refers to positive and negative terms. Only 
positive to negative transitions and vice versa are allowed. The s <=> s and a <=> a selection 
rules refer to symmetric and antisymmetric terms. Only symmetric to symmetric or 
antisymmetric to antisymmetric transitions are allowed. The g <=> u selection rule requires 
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only even states to odd and vice versa are allowed. These selection rules are subject to the 
limitations arising from compliance with coupling schemes to be covered in the next 
section (Hund's case (a) and (b)). 

TABLE 4.01 

APPROXIMATE MOLECULAR SELECTION RULES 
(Herzberg, 1950) 

AS = 0 
+ <=> - 
s <=> s 
a <=> a 
g <=> u 
AA = 0,± 1 
AJ = 0, ± 1 

There are exceptions to these selection rules. One exception of 
significance is that Z + to Ir is not allowed but Z + to Z + , and to Z' is allowed. Other 
exceptions will be dealt with on an individual basis in Chapter V. 

Violation of one of the selection rules in Table 4.01 does not mean 
the transition is rigidly forbidden. The transition may be allowed as magnetic dipole or 
electric quadrupole radiation which introduces additional selection rules. These rules are 
supplemental exceptions to the selection rules of Table 4.01. The Herzberg bands for 
molecular oxygen as well as the Vegard-Kaplan and Lyman-Birge-Hopfield bands for 
molecular nitrogen all violate at least one of the selection rules listed in Table 4.01 (see 
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Chapter V). Table 4.02 lists selection rules that apply for magnetic dipole and electric 
quadrupole radiation. 



TABLE 4.02 

SUPPLEMENTAL SELECTION RULES FOR ELECTRIC UUADRUPOLE 

A ND MAG N ETIC DIPQLE RADIATIO N 

(Herzberg, 1950) 

+ <=> + 

- - 

g«g 

U <£=> U 



B. SYNTHETIC SPECTRA 
1 . Overview 

The previous section defined molecular transitions along with the 
nomenclature and selection rules that apply to molecular transitions. In this section the 
wavelengths and relative intensities of these transitions will be identified. The synthetic 
spectra that are produced by these techniques can then be used to link observed emissions 
with the emitting molecule. Unique emission features allow observers to identify the 
"fingerprints" of each emitting species. The techniques described in this section can be 
used by computer modelers to develop synthetic spectra for each emission band that will be 
observed. 

2 . Line Positions 

The magnitude of the energy difference shown in equation 4.9, divided by 
he, for either an emission or absorption, is known as the line position. Line position 
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carries units of cm -1 . Another descriptive name is the wave number. The line position a , 
follows from the change in the appropriate parameters for the electronic, vibrational, and 
rotational contributions to the total. These are: the total electronic energy T e , the vibrational 
energy F(J), and the rotational energy G(V). The vibrational and rotational energies remain 
as were discussed earlier in Chapter IV. For a given transition from an upper state 
(T e ',v' J') to a lower state (T e ",v",J") , the line position a is expressed as is 

a = T - T" + G'(v) - G"(v) - F'(J) - F"(J). (4.10) 

The electronic portion of the entire band system a e = T - T" remains constant. This can be 
used to simplify Equation 4.10. The line position for a particular band system can be 
written as 



a = <y e + G’(v') - G"(v") - F'(J') - F"(J"). (4.11) 

For any v' to v" transition, there may be hundreds of rotational transitions. The wave 
number for any one of these lines becomes as is shown by:. 

Ofj'j" = a 0 + F'(J') - F"(J"). (4.12) 

The <7o term is known as the band origin or band head. Its value is given by 

<y 0 = <r e + G'(v') + G"(v"). (4.13) 



39 



At this point the rigid rotator and simple harmonic oscillator approximations 
need to be addressed. The constants used to determine line position are given by Huber 
and Herzberg (1979). However, these constants include anharmonicity and non-rigid 
rotator effects. In this case the vibrational term becomes 

G(v) = (o e ( v + j) - (D e x e (v + y). (4.14) 

The values for (Oe and x e are constants of each molecule and the second term accounts for 
anharmonicity. The rotational term becomes 

F(J) = B V J (J + 1) - D v J2(J + 1). (4.15) 

where the second term accounts for the non-rigid rotator. The factors B v and D v depend on 
the vibrational level. The expression defining B v is given by 

B v = B e - a e ( v + j) + higher order terms. (4.16) 



Similarly the expression for D v is 

D v = D e - Pe( v + ^-) + higher order terms. (4.17) 

The constants B e , a, D e , and P are all found in Huber and Herzberg. 

The spin Z * 0 and orbital angular momentum A * 0 contribute to molecular 
degeneracy through several coupling schemes. The rotational term values in these cases are 
highly dependent on the coupling between A and Z. In 1927, Hund (1927) first 
considered five possible cases of coupling between A and Z. The first two of these are 
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applicable to molecules with emission bands between 1800 and 3400 A and will be 
examined further. The total electronic angular momentum vector Q is the algebraic sum of 
the A and L components. The algebraic rather than vector sum shown in equation 4.18 is 
appropriate since both terms are aligned along the intemuclear axis. 

Q = | I + A | (4.18) 

Hund's Case (a) is appropriate when strong coupling exists between the 
electronic orbital and spin motion, and the intemuclear axis. At the same time electron 
motion is weakly coupled to the rotation of the molecule. This situation is completely 
analogous to a symmetric top. The electron motion around the two nuclei results in a non- 
zero moment of inertia about the intemuclear axis. The molecule will precess about the 
total angular momentum vector J, just as the symmetric top does. The symmetric top is 
shown in Figure 4.06. The curved arrow in the figure represents the rotation of the entire 
molecular system about the J vector. A indicates the orbital angular momentum of the 
electrons about the intemuclear axis. N is the component of the total angular momentum 
perpendicular to the intemuclear axis. 
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Vector Diagram for a symetric top 
Figure 4.06 The Symmetric Top (Her/.berg, 1950) 

The system in Hund's Case (a) has been complicated by the interaction of 
the spin and orbital electronic motion with the nuclear rotation and with the intemuclear 
line. The ground state, where spin and orbital angular momentum go to zero, will be 
of the same configuration as the symmetric top. Figure 4.07 illustrates the coupling in 
Hund's case (a). In case (a) the nuclear rotation and electronic motion coupling is weak, 
while the coupling between electronic motion and the intemuclear line is strong. The 
magnitude and direction of the J vector is constant, while the electronic angular momentum 
Q and the angular momentum of nuclear rotation, N, couple to produce the nutation about 
the total angular momentum J. Additionally L and S are precessing much faster than the 
precession about J. The bold arrow in Figure 4.07 indicates the nutation of the system 
about J while the ovals at the end of the A and X vectors indicate the more rapid precession 
of L and S about the intemuclear line. 
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Hund's Case (a) 

Figure 4.07 Hund's Case (a) (Herzberg, 1950) 

Hund's case (b) has a substantially different coupling scheme. The spin S 
is only weakly coupled to the intemuclear axis. The orbital part of the electronic motion is 
still strongly coupled to the axis but more moderately coupled to the rotational motion. 
Additionally, when A =» 0 case (a) can no longer apply. The K vector encompassing the 
total angular momentum apart from spin is then introduced. K is a quantum number with 
integral values as defined below: 

K = A, A + 1, A + 2, ... (4.19) 

The total angular momentum J now follows from 4.19 as 

J = K + S, K + S - 1, K + S - 2, . . . | K - S|. (4.20) 
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The system configuration and relative coupling for Hund's Case (b) are given in Figure 
4.08. 



N 



Hund's Case (b) 

Figure 4.08 Hund's Case (b) (Herzberg, 1950) 

The amount of splitting of the energy levels of a multiplet state depends on 
the molecule's adherence to Hund's case (a) or case (b). A general rule is that multiplet 
splitting is large for Hund's case (a) and approaches zero for Hund’s case (b). The amount 
of splitting can be accounted for by the expression 




T e = To + AAE . (4.21) 
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In equation 4.21 A and Z are as have already been defined and A is a molecular constant 
for the particular muliplet term. 

The Hund's coupling cases discussed above are idealized limiting cases. In 
general coupling does not neatly comply with one case or another. As rotational energy 
levels increase, a transition between case (a) and (b) is observed. This is known as 
uncoupling and there are two types of uncoupling mechanisms between case (a) and (b). 
These are 

A-type doubling and spin-type uncoupling. Both mechanisms will cause an energy level 
splitting. A-type doubling occurs for Hund's case (a) and (b), for A * 0 states (II, A, O 
etc.). This doubling is very small relative to the other energy differences involved in a 
transition. Figure 4.09 illustrates A-type doubling for two 2 I1 states. The splitting steadily 
increases as J increases. 




Figure 4.09 A-type Doubling for the 2 I1 State (Herzberg, 1950) 
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Spin-type doubling depends on the spin S for the molecule. As rotational speeds 
increase, the rotational speed becomes comparable to the precession of S about L. Hund's 
case (a) is dominant for small rotations and splitting is large. For large rotation the 
situation approaches Hund’s case (b) and splitting is reduced. The result of this larger 
splitting due to spin-type uncoupling results in different rotational term values for different 
spin sub-states. 

3. Molecular State Populations 

The calculation of the population of a specific molecular state requires 
accounting for all factors that contribute to the population. These factors vary between 
Hund's case (a) and (b). A state in Hund's case (a) is defined by n, v, and Z, while case 
(b) is defined by n, v, and N. For those states where A>0, A-type doubling introduces a 
parity quantum number p that must be accounted for. Although line positions are not 
sensitive to A-type doubling, the calculation of the the population rate and transition 
probability requires the two sub-states be considered. This is because each A-doubled 
component is individually populated. Transitions from each upper A-doubled state 
constitute a fully resolved line. 

The volume emission rate can be calculated by considering the fractional 
population of the upper state multiplied by a spontaneous emission coefficient between the 
two states. This is shown for Hund's case (b) in the equation below: 

E(n , v'N , J , p',n"v"N"J"p") = N(nVN'J'p') A(nVN'J'p\n"v"N"J"p"). 
(4.22) 

where E is the volume emission, N is the fractional population, and A is the Einstein 
coefficient for spontaneous emission between two states. The factors enclosed are 



46 



parameters defining the electronic state n, the vibrational level v, the rotational level N, the 
vector J, and the parity p. For Hund's case (a) the parameter I would replace the rotational 

level N. 

The number of molecules in a specific n' and v' level may be accounted for 
by the production rate divided by the loss rate in the steady state using 



N(n'v') = . (4.23) 

^ Ay'y M 

Now Equation 4.22 can be multiplied by 1 in the form of N(n'v')/N(nV) and expressed as 



E(n VN'J'p',n" v"N" J"p") = P(nV) N( ^^ ) ' P ' ) A J ’ J " . 

^A v ' v " 



(4.24) 



In equation 4.24 the factor Aj-j" is equal to A(n VN'J'p',n" v"N" J"p"). 

The production mechanism leading to emission for allowed transitions is 
photoexcitation. The production rate P(n V), is proportional to the solar flux and 
effective cross section for the transition to a specific electronic state. It is assumed that the 
upper state rotational population is statistically similar to that of the ground state. The 
branching ratio into various vibrational levels are given by the Frank-Condon factors. 

These factors arise from the Frank-Condon principle. A brief description of the principle 
follows. The line strength for a particular transition is a function of the relative positions of 
the wave functions between the upper and lower states. The electronic transition time for a 
molecule is on the order of 10~ 16 seconds. Because the vibrational transition time of 10' 14 
seconds is relatively long the intemuclear distance will remain virtually constant during the 
electronic transition. Referring back to Figure 4.04, this transition may be expressed as a 
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vertical line between electronic states. The maximum probability density in the excited state 
will determine the upper state location. A strong wave function spatial overlap will 
determine the most likely lower eigenstate. The most favored location is generally at the 
end of an oscillation where the eigenfunction is a maximum. The relative strength of a 
transition from a vibrational level (v') of an excited electronic state to the vibrational level 
(v") of a lower electronic state is then based on the locations of the upper and lower levels, 
and the probability amplitude at the respective upper and lower locations. This concept, is 
known as the Frank-Condon principle (Eisberg, 1985). 

The production rate for an allowed transition can be written as 

np2 

P(n\v’) = n fxn' q„v- (4.25) 



where rtF is the solar flux (photons cm* 2 sec' 1 A* 1 ), X is the wavelength of incoming 
radiation, f*ri is the oscillator strength for the electronic transition from the ground state to 
the n' state, and q ov ' is the Frank-Condon factor between the v=0 level of the ground state 
and the v' level of the upper n' state. This production equation makes the valid assumption 
that virtually all molecules are in the ground state prior to absorbing incident radiation. 



Tatum (1967) develops a detailed description to determine the relative 

populations of molecular energy levels. The development is summarized here. The 

Nln'v'N'J'p') 

fractional population N ( n ' y ' ) ^ ° * a ®* ven state 8 enerall y obeys a Boltzman 

distribution. The Boltzman distribution for a state can be expressed as 



Ni 

N 



Xgje-Ej/kT 

j 



(4.26) 
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In Equation 4.26, gi is the statistical weight of the ith group and the denominator is the 
associated partition function. Tatum applies this distribution to the different molecular sub- 
states corresponding to appropriate Hund's coupling cases. As an example for Hund's 
case (b) in states where A * 0, the fractional population may be expressed as 



N(n'v , N'J'p') he B v ^ (2J + 1) 
N(n’v') " kT ®(2S + 1) 



exp( 



hcF(N) 
kT > 



(4.27) 



The factor <I> introduced in Equation 4.27 is a constant for hetero- or homo- nuclear 
molecules. The various other factors in equation 4.27 have all been defined previously. A 
similar expression is obtained for Hund's case (b). 

A i ? i** 

The remaining factor in equation 4.24 — ■ LJ — , is yet to be determined. 

The factor Ajj" is the transition probability. This probability for the rotational branching 
between two vibrational levels can be determined by the Honl-London factors. These 
factors sum to the value 2J +1 so in order to normalize the factor it must be divided by 2J 
+1. This will yield the correct branching probability. The Honl-London factors yield the 
probability of a rotational branch for a given v'-v” transition. To consider a transition on 
the specific v" level another factor A v 'v” must be included. The factor Aj'j" becomes 

Si!' I" 1 

Aj'j" = A vV ’ 2 ~ p ~ T - l * ( 4 - 28 > 




Equation 4.24 can now be rewritten as 



E(nVN'J'p',n"v"N"J"p") = 



N(nVN'J’p') 

N(n'v') 



A v 'v" 




S(J’J") 

2 J ' + 1 * 



(4.29) 
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The single-scatter albedo ay v " of a fluorescent band is defined as 



CDv'v" 



A v * V ** 




(4.30) 



The factor A v V" is the transition probability of a molecular band which is related to the band 
strength by the relation shown below: 



A v'v" 



64jt4v 3 v 'v" 

3h~c 3 dT 7 



R 0 v'v". 



(4.31) 



In Equation 4.31 V v 'v" is the frequency of the band, d v ' is the degeneracy of the upper 
level and Re 2 qv'v" is a close approximation of the band strength S v v- These factors 
allow Equation 4.28 to be expressed as 



<3v’v" 



V 3 v’v 



n 




v"qv'v" 



(4.32) 



Equations 4.30 through 4.32 can now be substituted into Equation 4.29 to yield the volume 
emission rate. 

Production for forbidden transitions conforms to the required calculations 
for electron impact excitation. The relative speed and abundance of electrons in the 
ionosphere results in collisions between electrons and neutral molecules being the dominant 
production source for forbidden transitions. Excitation by electron impact becomes the 
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dominant source of the excited state for the transition. The source equation is now highly 
dependent on the electron flux, as well as the electrons cross section for excitation. The 
emission rate as compiled by Barth (1965) is given by 

Emax 

4^ / v 'v" = J^Ie^e dE q ov ’ (OvV No ch(9) (4.33) 



where Ie is the particle flux in electrons cm' 2 sec* 1 . Both, the electron particle flux and the 
electron cross section for excitation of an electronic state <Je are energy dependent. The 
Franck-Condon factor for excitation q ov >, single-scattering albedo cbvV', and molecular 
column density A j q, adjusted by the Chapman function ch(0), remain consistent with 
those for the photoexcitation case. It is necessary to integrate over all available electron 
energies when a range of electron energies is present. This is the usual case in the 
ionosphere. 

The methods of this section have been used to compile the synthetic spectra 
to be illustrated in Chapter V. Each spectrum is highly dependent on the unique properties 
of each molecule involved. Separate algorithms have been used for each emission feature 
in Chapter V. The actual calculation requires that a volume emission rate for each line in 
the band be calculated. The contribution from each bin is then summed to produce the 
spectrum. An example of this process is shown in Figure 4.10. In this figure the (1,0) 
gamma band of Nitric Oxide is shown for a wavelength bin size of .03 A. The individual 
lines are distinct. The spectra in Chapter V have a lower resolution comparable to 
instrument resolving capability. As a result individual rotational lines are not 
distinguishable. 
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Wavelength (X) 



frigure 4.10 Synthetic Spectrum for (1,0) Gamma Band Emissions 

(Cleary, 1985) 
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V. DAY GLOW EMISSION FEATURES IN THE 1800-3400 A 



WAVELENGTH RANGE 



This chapter is concerned with identifying observable emission features in the 
wavelength range of the instrument. N 2 , 02, O and NO all have emissions in the 1800- 
3400 A window. The transitions result from the interactions presented in Chapter IV. 

There are additional emission bands that occur in the window besides those that will be 
described. These are of such a low intensity at the altitudes of interest that they are not 
measurable. The bands and lines presented in this chapter are primarily the ones that are 
predicted to be observed, during the flight. Computer simulations of synthetic spectra have 
been created by McCoy (1981), Cleary (1985), and Conway (1982). 

The observed intensity will be limited by the detector efficiency. The instrument 
detector efficiency is limited at the long end of the experimental wavelength range by 
photocathode efficiency. Figure 5.01 shows this efficiency over the wavelength range of 
the instrument. The solid curve in the figure demonstrates that detector efficiency drops off 
rapidly past 3000A . This efficiency drop will become a significant consideration when 
attempting to identify emission features of longer than 3000A that are observable. 
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Photocathode 

Efficiency 




Figure 5.01 Photocathode Efficiency 



A. MOLECULAR NITROGEN 

Molecular nitrogen has three observable band systems: the Vegard-Kaplan (A^Z U + 
- XlZ g +), 2nd positive (C^riu - B^FIg), and the Lyman-Birge-Hopfield (a^IIg - X^rig) 
bands. The Birge-Hopfield (b^n u - X^Z+g) band was also investigated as a possible 
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contributor but its intensity in the experimental window was found to be negligible. The 
potential energy curves for N 2 are shown in Figure 5.02. References to specific energy 
levels of the N 2 molecule will be from this figure. In addition, Figure 5.03 is a schematic 
energy level diagram of N 2 . This figure illustrates the major transitions of molecular 
nitrogen. It also displays the relative energy differences for the emission bands. 

It is of note that molecular nitrogen is the dominant component of the atmosphere at 
the lower end of the experimental altitude range. However, at around 175 km atomic 
oxygen becomes the dominant species (see Figure 2.08). Since N 2 is the dominant 
species, a reaction involving N 2 may be observable although it has a relatively low 
emission probability. Additionally, the relative importance of N 2 is strongest at lower 
altitudes and will tend to drop as altitude increases unless some other factor arises such as 
quenching at lower altitudes. The N 2 emissions in the ultraviolet tend to be weak because 
they are a result of forbidden transitions subject to quenching, or a combination of both 
factors. 
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POTENTIAL ENERGY (dec Iron vo. 




Figure 5.02 N 2 Potential Energy Curve 
(Journal of Physical Chemistry Reference Data, 1977) 
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Figure 5.03 N 2 Energy Level Diagram (Sharp, 1970) 

1 • Vegard-Kaplan System 

The V- K system represents the A3l u + - xll g + transition, and has a range 
of 1250 to 5325A (Krupenie, 1977). In 1932, Vegard first discovered the system by 
bombarding solid nitrogen with electrons. Later in 1934, Kaplan succeeded in obtaining 



57 



V-K bands with gaseous nitrogen. Kaplan also adjusted the wavelengths of the lines 
proposed by Vegard. (Krupenie, 1977) 

Barth (1965) identified both fluorescence and electron processes as 
production sources for the V-K emissions. However, Sharp (1970) has determined that 
the sources for the upper state are almost exclusively due to electron impact. The upper 
state A 3 L + U is populated as a result of direct excitation by electron impact, as well as 
cascading from higher states (Broadfoot and Hunten,1964). Dominant sources of cascade 
are directly from the B^flg state and indirectly from the C3fl u state. It is of note that the 

relative population of the vibrational levels depends on the excitation source. Table 5.01 
has been computed to demonstrate the relative distribution of vibrational levels as a result of 
each case. The first row considers the distribution for direct excitation to the A 3 Z + „ state. 
The second row considers the relative distribution of vibrational states only for cascade 
from higher states. While direct excitation results in a concentration at higher vibrational 
levels, cascade results in a concentration at the lower vibrational levels. Table 5.01 is in 
non-normalized units where a larger number indicates a higher relative population of the 
entire A 3 Z + „ state. 



TABLE 5.01 

RELATIVE POPULATIONS OF THE A3£+„ _STATE 



(Ahmed, 1969) 
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This scheme is further complicated by the proposal by Campbell and Thrush (1967) and 
later by Gilmore (1969) that inverse first positive transitions act as non-radiative loss 
mechanisms (see Figure 5.03). It has been suggested by Torr and Torr (1982) that the 
complex population mechanisms for the V-K system require further data. This experiment 
may be able to provide some useful data at the upper end of the experimental altitude range. 

Because the upper state of the Vegard-Kaplan state is metastable, quenching 
rather than emission will dominate. This is particularly true at the low altitude end of the 
experiment. It has been shown that the Vegard-Kaplan band is a weak feature in the 
dayglow although it is quite bright in the aurora (Torr and Torr, 1982). The primary 
reason for this appears to be quenching. Shemansky (1969) determined that the upper state 
possessed a mean lifetime of 1.8 seconds. The atmospheric collision frequency is about 
2.25 per second at 120 km ( CIRA model 1965, in Sharp, 1970). Therefore at low 
altitudes there is a likelihood of four collisions prior to an emission. Quenching may result 
from collisions between N 2 and neutral atoms or molecules. Candidates include N 2 , O 2 , 

O, NO, and N as potential quenching agents (Torr and Torr, 1982). Sharp (1970) has 
identified O as the dominant agent. The number density for all species drops by two orders 
of magnitude between 100 and 200 km (see Figure 2.08). Collision frequency can be 
expected to drop off proportionally. Thus, the Vegard-Kaplan bands will assume an 
increasingly significant role as altitude increases. 

Conway (1985) has identified the (0,9) band at 3353A as having an 
intensity of 2 kR at an altitude of 1 50 km. This is a detectable amount although it is well 
past the 3000A threshold of reduced detector efficiency. It is necessary to investigate other 
vibrational bands in the V-K system to determine whether or not this band system will be 
observable. There is a convenient factor to distinguish the relative intensities of emission 
bands for a particular band system. This factor is a combination of the Frank-Condon 
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factor for excitation and the albedo for scattering (Q- omega). This is interpreted as that 
portion of the total excitation to a particular electronic state that is excited to a specific v' 
state and emitted back down to a particular v" state. The Q.omega factor for the (0,9) 
feature is 2.47xl0' 5 . Numerous emission features below the 3000A cutoff have similar or 
even higher Q-omega factors (Barth, 1965). The highest is the (5,1) band at 184lA, with 
Q-omega equal to 1.48x10"^. This would support the contention that N 2 V-K should be 
observable in measurable quantities during this experiment. 

The V-K transition is spin forbidden but remains electric dipole radiation 
(Conway, 1980). It violates the As = 0 requirement as the spin 3 to spin 1 states, are not 
allowed. However,the 3£ u + to 1 Z g + transition will occur despite it being forbidden. 
Herzberg (1950) states that the spin orbit interaction causes the state to "borrow" some 
of the characteristics of the allowed 3 rio state, or the state "borrows" some of the 
characteristics of the *11 state. Four branches will result. The synthetic Vegard-Kaplan 
spectrum between 1800-3400 A is shown in Figure 5.04. 
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Wavelength (angstroms) 

Figure 5.04 N 2 Vegard-Kaplan Bands 

2. N 2 Second Positive System 

The second positive system represents the C^n u - B^rig transition. It 
ranges from 2680 to 5460A (Krupenie, 1977). The Second Positive system was first 
observed in 1869 (Krupenie, 1977). It appears in most nitrogen sources including pure 
nitrogen and air. It was observed in the day glow as early as 1945 (Pearse and Gay don, 
1963). Numerous observations include a recent experiment by Conway (1980) in 1980. 

Conway and Christensen (1985) reported an observed intensity of about 
lOkR for the (0,0) feature at 3370A during solar maximum conditions. It appears that this 
is the brightest feature in the experimental wavelength region. Other features in the region 
are less intense. The Q-omega factor for the (0,0) band is .280 . The highest factor below 



61 



3000A is the (2,0) with a value of .021 (Barth, 1965). This is an order of magnitude lower 
than the (0,0) band and is an indication that the 2nd-Positive band will make little 
contribution to the observed emissions in the experimental wavelength region. 

The intensity can be calculated using the following method. From Equation 
4.33, the intensity due to electron impact excitation is equal to the product of the total 
electron impact cross-section, slant column density, and Q-omega factors for the band. 

The calculation used the electron impact cross-sections for the C^lly state which is shown 
in Figure 5.05. The cross-section at each energy can be multiplied by the electron flux at 
that energy. The electron flux has been calculated in a simulation by Strickland and Meier 
(1982) using an F10.7 equal to 206, a solar zenith angle of 56 degrees, and a temperature 
at infinity of 1000K. The data points on the figure indicate separate results by E. Oran 
(Oran, 1978). The electron flux at each energy can be multiplied by the cross-section at 
each energy the result is shown in Figure 5.07. The integration over the electron energy 
range shown in Figure 5.07 will yield the total electron excitation cross-section for each 
state. This cross-section multiplied by the Q-omega factor and a representative column 
density at a given altitude will give the expected observable intensity in Rayleighs. The 
slant column density is computed by adjusting the vertical column density by the chapman 
function for looking across the limb. This has been done for the values given in Table 
5.02. 
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TABLE 5.02 



INTENSITY FACTORS FOR N?(C3lL. - B3n g )I.RAMSlima 

Em ax 

Ja E Ie dE = 8. 25xl0- 9 

% v '<jt>v'v" (Q-omega) = .28 (Barth, 1965) 

N =1.8xl0 18 mol crrr 2 { for N 2 n( 150km) = 1.0xl0 10 mol cm'3,ti n f= 1800°K, 

ch(30,90°)) 

AnI = 4.2kR 

This calculated intensity is a reasonable approximation to the measured value of lOkR, 
considering that both N 2 column densities and electron fluxes are approximated. 
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Electron Impact Cross-Section (cm2) 




Figure 5.05 Electron Impact Cross-Section vs. Electron energy N 2 (C 3 n„) 
and 02(A 3 Z + u ) states (Cartwright et al, 1977 and Borst, 1972). 
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Figure 5.06 Electron Energy vs Electron Flux (Strickland & Meier, 1982) 
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eE vs.(<Je * Ie) * for N2(C 3 n„) and C>2(A 3 I + U ) 




Figure 5.07 Electron Energy vs. (Cross-section x Electron Intensity) 

Because of the the AZ = 0 selection rule discussed in Chapter V, three sub- 
bands will result for Hund's case (a) 

3 n 0 - 3 n 0 
3 ni -3m 
3 ri2 - 3 ri2 

(Herzberg, 1950). Each sub-band has a strong P and R branch. The Q branch is a weak 
branch except for the 3 rio - 3 Ilo case. This results in triple band heads which are degraded 
to the blue. The Physical Chemistry Reference Data Manual (Krupenie, 1977) has 
identified fifty of the bands. When J is small, the transition conforms to Hund’s case (a) 
(see Figure 4.07). As J becomes larger, the coupling approaches that of Hund's case (b) 
(see Figure 4.08). This transition exhibits some A-type doubling doubling, as required by 
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Hund coupling for a transition. The second positive synthetic spectrum that will 

contribute to the experimental window is shown in Figure 5.08. 




Wavelength (Angstroms) 

Figure 5.08 N 2 Second Positive Bands 
3. Lyman-Birge-Hopfield Band 

The Lyman-Birge-Hopfield bands (a^g - X^+g) have a wavelength 
range from 1 160 to 3020 A (Barth, 1965). The emission bands were first observed in 
1928 by Birge and Hopfield (Pearse and Gaydon, 1963). They are readily visible as both 
emission and absorption features in the vacum ultraviolet. 

The transition is forbidden as electric dipole radiation Emissions will result 
from electric quadrupole and magnetic dipole transitions. Furthermore, the upper state is 
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forbidden from the ground state (Herzberg, 1950). The excitation will occur primarily as a 
result of electron bombardment. 

LBH emissions have been used as an accurate measure of the photoelectron 
flux (Conway, 1982 ). Conway observed dayglow intensities approaching 800 R at 160 
km for the (2,0) feature at 1383 A. The Q-omega factor for the (2,0) band is 4.63xl0' 2 . 
The (2,8) band at 1837A has a Q-omega factor of 1.87xl0" 2 . The (3,9) at 1854A, (3,10) 
at 1928A, (4,10) at 1870 and (4,11) at 1945A emission bands all have similar Q-omega 
factors; on the order of 10‘ 2 [see Barth (1965) for the complete table electron impact cross- 
sections]. This would indicate that LBH emissions should be visible in measurable 
intensities in the wavelength region. These features are at the short wavelength end of the 
detector which has a relatively low efficiency . Nevertheless, it would seem that the 
Lyman-Birge-Hopfield Bands should contribute to the observed experimental intensity. 

The peak observed intensity for all LBH emission features should be in the range of 150 to 
200 km in altitude if atmospheric conditions during Conway's observations remain 
generally in effect. 

The a^g - X^+g transition violates the g <=> u selection rule given in Table 
5.1. The transition is allowed by the special selection rules of Table 5.2. The AJ selection 
rule for quadrupole radiation becomes 

AJ = 0, ±1, ± 2 (Herzberg, 1950). 

The resultant branch structure is some what complicated. It is best presented graphically as 
is shown in Figure 5.09. 
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Figure 5.09 n g - 1 Z + g Branches for Magnetic Dipole and Quadrupole 

Bands (Herzberg, 1950) 

The branch structure of the LBH system will produce the lines shown in Figure 5. 10 in the 
experimental window. 




WavElength (Angstroms) 

Figure 5.10 N 2 LBH Bands (Cleary, 1985) 



69 



MOLECULAR OXYGEN 



B . 

There are two observable bands of O 2 that may contribute to the experimental 
window. These are the Herzberg I, and Schumann-Runge (SR) bands. The first of these 
is forbidden, while the second is an allowed transition. The potential energy level diagram 
for O 2 is shown in Figure 5.1 1 . The schematic energy level diagram for O 2 is shown in 
Figure 5.12. This figure demonstrates the major transitions of molecular Oxygen. It also 
displays the relative energy differences for the emission bands.The Schumann-Runge 
bands have an easily accessible threshold energy for photo dissociation which minimizes 
the observed intensity. The energy curve for the B-^Z'u is of note. This is the upper state 
for the Schumann-Runge transition. O 2 in this state will dissociate into O(-P) + OC 1 D) 
atoms. It may also predisscociate to the O(^P) + O(^P) atoms. The upper state for the 
Herzberg system (A- Z + u ) may also result in the O(^P) + 0(3p) atoms. 




Figure 5.11 O 2 Potential Energy Diagram (from Gilmore, 1965) 
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Figure 5.12 O 2 Energy Level Diagram 
1 . Sehumann-Runge Bands 

The Schumann-Runge system (B 3 Z* U - X 3 l* g ) ranges from 1750 to 5350 
A. With over 100 single-headed red degraded bands, the Schumann-Runge (SR) is a very 
extensively studied Oxygen system. It was identified by Schumann as early as 1903 
(JPCRD, 1972 ).The large number of lines arises from the greater intemuclear distance of 
the upper electronic state over the lower state. This results in a large number of accessible 
lower vibrational levels for each upper vibrational level. 

The B 3 Z' U state has a strong absorption cross section However as 
mentioned above, this state lies very close to the Oxygen dissociation continuum. The 
ionization and dissociation continuum are shown with respect to the absorption cross 
section in Figure 5.13. 
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Figure 5.13 Absorption Cross-Section for O 2 , Showing 
Continua, and SR Bands (Banks, and Kockarts, 1973) 

For v' levels greater than three, the Cb molecule will predissociate rather 
than fluoresce (Banks and Kockart, 1973 ). This means only the v'=0,l,2,3 states will 
contribute to the emission spectrum. In fact the Schumann-Runge system is actually the 
major source of atomic oxygen. 

An estimate of the SR intensity can be made by comparing its g factor and 
the O 2 slant column density with the g factor and slant column density of known 
emissions. Barth (1965) lists the number of photons scattered per second per atom (g 
factors) for numerous atmospheric constituents. These g factors multiplied by the slant 
column density determine the relative strength of an allowed transition. The g factors for 
the first four v' states for the Schumann-Runge bands are extremely small. The largest 
value in the experimental wavelength range is 3.59X10' 11 for the (3,7) band at 2442 A. To 
put this in perspective a bright allowed feature such as the Nitric Oxide (1,0) gamma band, 
has a g factor of 7.7x1 0- 6 . This implies that an NO Gamma transition is much more likely 



72 



than an O 2 Schumann-Runge transition. However, a contrasting value from Figure 2.08 
shows that the number density for O 2 is about five orders of magnitude higher than NO. 
This would seem to imply that the strongest O 2 SR bands should be of comparable 
intensity to NO gamma bands. The values used to calculate the (3,7) SR band at an altitude 
of 150km for active solar conditions is provided in Table 5.03. 

TABLE 5.03 

ll2_SRn-7tBAND INTENSITY CALCULATIONS 

g = 8.25x10-9 

N =1. 8x10*8 mol cnr 2 {for O2n(150km) = l.OxlO^mol cm -3 , tinf = 1800°K, 

ch(30,90°)} 

4 n I = 4.2kR 

The Oxygen Schumann-Runge band system was not identified by Sharp (1986) in the 
thermosphere between 2000 and 8000 A. No other experimental observation of this band 
has been found. However, it appears that emissions do exist in the experimental 
wavelength region, and they may contribute to the observed intensity. 

The 3 Z‘ U state is an allowed state for photon excitation from the X 3 Z~g 
ground state. Once in the excited state the 3 Z' U - to 3 L‘g transition produces both P and R 

branches. Each branch has three components for a total of six main branches for each 
band. Additionally there are six satellite branches where AJ * AK. Splitting will become 

significant only when K becomes large (Herzberg, 1950 ). The O 2 molecule has zero 
nuclear spin. This condition requires the molecule to follow Bose statistics (Herzberg, 
1950 ). This requirement produces a band structure where every alternate triplet in the P 
and R branches appears. The band system is red degraded. 
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The synthetic spectrum for the SR bands that will fall in the window is 
shown in Figure 5.14. The synthetic spectrum shown is for all the v'=0, 1,2, 3. For v'>3 
the molecule will predissociate rather than fluoresce. 




Wavelength (angstroms) 



Figure 5.14 O 2 Schumann-Runge Bands 
2 . Herzberg Bands 

The Herzberg transition (A 3 X + U - X 3 Z' g ) appears to range from 2590 to 
4880 A (JPCRD, 1972 ). The band system appears to have been discovered rather recently 
by Herzberg, in 1952 (Pearse and Gaydon, 1963). The band system is a prominent feature 
in the nightglow, where the excitation source is inverse predissociation (Barth, 1965 ). 
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As a forbidden transition, Herzberg bands will result primarily from 



photoelectron impact. Its dayglow intensity will be less than a corresponding transition in 
the Schumann-Runge Band. The Z + <=> Z + or L* <=>Z • selection rule is the selection rule 
for this transition. 

The upper state is subject to dissociation for incident energies greater than 
5.1eV. It corresponds to a wavelength of 2420A. This means that at any altitude, 
emissions as a result of the A 3 Z + U - X 3 Z* g transition will be longward of 2420 A. The 
Herzberg dissociation peaks at 40 km and is closely linked with ozone production 
(Chamberlain, 1987). The A 3 Z + U state dissociation reaction is shown below: 



The calculation to determine the intensity of the Herzberg bands follows 
from that developed for the N 2 Second Positive bands. Figures 5.05, 5.06, and 5.07 are 
used to determine the intensity. The calculations were done for active solar conditions, 
looking across the limb. The brightest band identified is the (12,1) band at 2528A. The 
factors used to determine the observed intensity are shown in Table 5.04. 



N =1.8xl0 17 mol cm* 2 {for 0 2 n(150km) = l.OxlO 9 mol cm* 3 , 1^= 1800°K, 



0 2 + e*(E > 5.1eV) => 0( 3 P) + 0( 3 P) + e. 



TABLE 5.04 



INTENSITY FACTORS FOR Q z (Herzbeni)TRANSITION 




%v'GVv" (Q-omega) = 1.98xl0* 3 (Barth, 1965) 



ch(30,90°)} 



4n/ = 45R 
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The calculated intensity of 45R is a marginal value as an observable value in the wavelength 
region. The Herzberg bands will contribute to the total observed intensity, however 
separate emission features may not be distinguishable. The O 2 Herzberg spectrum that 
contributes to the experimental window is shown in Figure 5.15. 




Wavelength (Ang stroma) 

Figure 5.15 O 2 Herzberg Bands 

C. NITRIC OXIDE 

Nitric Oxide has three band systems in the experimental window. The Gamma 
band (A 2 I + - X 2 I1), Delta band (C 2 n - X 2 II), and Epsilon bands (D 2 I+- X 2 n>. The 
potential energy diagram for NO is shown in Figure 5.16. A schematic energy level 
diagram, which illustrates the major NO transitions, is shown in Figure 5.17. 
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Although NO has a significantly lower density compared with N 2 , O 2 , and O, the NO 
bands are very intense. As shown in Figure 2.08, at the lower end of the experimental 
altitude it is approximately five orders of magnitude less than the N 2 and O 2 densities. 
However, all of the the NO band systems are allowed transitions. The importance of nitric 
oxide as an emission source of the airglow becomes as significant as a corresponding 
forbidden transition of a more abundant species. The importance of NO diminishes with 
altitude beyond its peak number density at about 105 km. A fourth system, the 
Beta band (B 2 n - X 2 n), also has transitions in the experimental window. However, the 
spatial overlap for this system results in intensities four to five orders of magnitude less 
than the other bands. It is not a significant contributor to the window. 




0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.2 3.6 



INTERNUCLEAR DISTANCE^) 

Figure 5.16 NO Potential Energy Diagram (Cleary, 1985) 
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Figure 5.17 NO Energy Level 



Diagram 
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1 . Gamma Band 

The Gamma Band (A 2 Z + - X 2 !!) is the most intense of the NO transitions 
in the experimental wavelength range. The bands range from 1873 to 6126 A (Barth, 
1965). The band system was observed as early as 1917 by Strutt (Pearse and Gaydon, 
1963). It commonly occurs in flames containing nitric oxide as well as shock-treated air. 
The (1 ,0) Gamma band, at 2149 A, is the most intense feature of Nitric Oxide in the 
ionosphere. The Gamma bands were not observed in the ionosphere until 1964 by Barth 
(1965). 

The primary excitation mechanism for the upper state is photoexcitation. 
Other sources to the A 2 Z + state include cascade from the C 2 Z + state(Sharp, 1986) and the 
chemical reaction shown below: 

NO(B 2 ri; v=7) + M =» NO(A 2 I+; v=3) + M(Young & Sharpless, 1962) 

These sources result in emission, where the NO molecule returns to the X 2 !! ground state . 
The upper state has a lifetime of 2x1 O' 7 seconds (Callear and Pillings, 1970) which insures 
emissions rather than collisions will dominate activities in the state. 

Measurements of the most intense Gamma bands are limited by self- 
absorption. Self-absorption has been observed by Cleary (1985) below 125 km for the 
(1,0) band at 2 149 A, and below 115 km for the (2,0) band at 2047A. These two bands are 
the strongest emission features of the gamma band system. Self-absorption is proportional 
to the column density times the oscillator strength for a band. This is why self-absorption 
is most significant for the strongest bands. 

Nitric Oxide Gamma bands have been observed recently by Cleary (1986) 
in 1982, during active solar conditions. His observations yielded observed intensities of 
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about 30kR for the (1,0) band at an altitude of 105 km. This intensity had diminished to 
2kR by an altitude of 190 km. Another measurement of the NO (1,0) gamma band at 2150 
A by Thomas (1978) was conducted during solar minimum conditions. His measurements 
indicate a measured slant intensity of 10 kR at 100 km, diminishing to less than 100R at 
300 km. These figures support the prediction that the NO Gamma band will be a 
significant contributor to the experimental window, especially at low altitudes. 

The upper state A 2 Z + , conforms to Hund's case (b) while the ground state 
is an intermediate case of (a) and (b). The 2 Z - 2 n transition results in doublet splitting for 
both the upper and lower states. This produces a total of four band heads. Additionally, 
the AJ =0, ± 1 selection rule will result in P,Q, and R branches. Therefore the entire band 
will contain twelve branches. Bands are degraded to the blue (Herzberg, 1950). The 
emission spectrum of the NO Gamma band emission spectra is shown in Figure 5.18. 
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NO Gamma Band 




Wavelength (Angstroms) 



Figure 5.18 NO Gamma Bands (Cleary, 1985) 

2 . Epsilon Bands 

The NO Epsilon bands (D 2 Z + - X 2 I1) are the same type of electronic 
transition, 2 Z + - 2 n, as the gamma bands. The band system ranges appears to be most 
intense in the ultra-violet although precise wavelength ranges have not been identified. 
This emission system is one of the lesser known in the NO family, although it was 
observed by Guillery in 1927 (Pearse and Gaydon, 1963). Separate emission features in 
the ionosphere were identified by Cleary (1986). 

The emission is a similar resonance scattering as the gamma bands. This means 
that an excitation from the ground state to the excited state leads to a relaxation back to the 
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ground state. The exclusive excitation mechanism for the D 2 Z + state is photoexcitation. 
(McCoy, 1981 ). 

Cleary (1986) measured a peak intensity approaching 700 R/A for the (0, 1) band at 
1950 A at an altitude of 105 km. It appears that this band as well as the slightly weaker 
(0,2) at 2020A and (0,3) at 2090A features should be observable in the experiment. 

Once again each band has four band heads, which yields the same twelve branch 
structure as the gamma bands. As shown in the potential energy diagram, an epsilon 
transition will be generally of higher energy than a corresponding gamma transition. The 
NO epsilon bands that will contribute to the experimental wavelength region are shown in 
Figure 5.19. 




Wavelength (Angstroms) 

Figure 5.19 NO Epsilon Bands (Cleary, 1985) 
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3 . NO Delta Bands 

The Delta Bands (C 2 n - X 2 II) are also an allowed transition for the 
2 n - 2 n case. Delta bands range from 1621 to 4288 A. Early observations of this band 
system were made by Knauss in 1928, and Schmid in 1930 (Pearse and Gaydon, 1963). 

The Delta excitation mechanism is photoexcitation similar to the Epsilon and 
Gamma bands. Delta bands do not demonstrate a fluorescence efficiency of unity, 
however. The accepted efficiency factor is .25 (Cleary 1986 ). This low efficiency is 
required because of predissociation. 

Cleary (1986) also observed delta bands during his 1982 experiment . 

Sharp (1986) has not identified any more recent observations. The strongest feature 
observed by Cleary was the (0,1) band at 1980 A. He measured a peak intensity of about 
150 R/A which will be difficult to distinguish in the experimental wavelength region. 

Once again each band has four band heads, which yields the same twelve 
branch structure as the gamma bands. As shown in the potential energy diagram, a delta 
transition will be generally of higher energy than a corresponding gamma transition, 
although it is of a slightly lower energy than the epsilon bands. The NO delta bands that 
will contribute to the experimental wavelength region are shown in Figure 5.20. 
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Wavelength (Angstroms) 

Figure 5.20 NO Delta Bands (Cleary, 1985) 

D. ATOMIC OXYGEN 

There is one atomic oxygen line that occurs in the experimental window. This is 
the 2972 A line. The line is a result of the (0(*S) - 0( 3 P)) transition. It is a well known 
forbidden emission of the oxygen atom. 

This line is represented by the transition shown below. 

0(!S) => 0( 3 P) + lvu(2972A) 
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This transition is forbidden as it violates AS = 0 requirement for atomic transitions. 
Although forbidden, this transition has been reported by Sharp at an average intensity of 
80R. This is a marginal observed intensity. Furthermore, this line is near the 3000A 
instrument efficiency cutoff which will tend to further limit the observed intensity. 
However, the relative density rises steadily throughout the experimental altitude range. The 
2912k line will be swamped by stronger emissions at low altitudes in the experimental 
region but could be distinguishable at higher altitudes. This will be especially true when O 
becomes the dominant constituent above 175 km. 

The upper state for this transition, the O(^S) state, is the source of two branches. 
These two branches are shown in the energy level diagram of Figure 5.21. 



SINGLETS TRIPLETS QUINTETS 




Figure 5.21 O Energy Level Diagram 
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It is clear from Figure 5.21 that an O atom in the l S state can either emit a 2972 or 5577 A 
photon. Furthermore, since the branching ratio between the 5577 and 2972 A lines has 
been determined a measurement of the 2912k. feature can be used to determine source terms 
for the O^S) state. The sources for the 0( 1 S) state are shown in Table 5.05. 

TABLE 5.05 

QXX£EKJ}S1 SOURCES 

O + e 
0 2 + + e* 

o 2 + + n 

N 2 (A3z u +) + O 

These sources include N 2 in the upper state for Vegard-Kaplan bands. This relationship 
could be used to assist in determining the quenching rate for the A 3 Z + U state of N 2 . Cross- 
sections and efficiencies for these source terms are given by Sharps (1986). 

The 2972 A line should be observable in the experimental window. It is near the 
limit of optimum detector efficiency, but below the 3000 A threshold of reduced efficiency. 
It appears that this line may also be useful in determining N 2 molecular densities as well. 
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VI. CONCLUSION 



A. SUMMARY OF FINDINGS 

Identification of emissions in the wavelength range of the MUSTANG experiment 
was accomplished using experimental results to support calculated emission feature 
intensity. All emission bands that will contribute to the experimental spectrum, as well as 
those features that will be distinguishable, have been tentatively identified. These results 
confirm the contention that the instrument wavelength range allows observation of all 
midultraviolet airglow emission bands. The extended wavelength and altitude range for 
this experiment allows a unique experimental opportunity to observe all major neutral 
atmospheric constituents in one spectrum. 

The major midultraviolet emission bands are shown in Table 6.01 . This table also 
identifies the brightest emission feature, its wavelength, and whether or not the band will 
be distinguishable. The last column provides special circumstances for observations of that 
band system. This column considers the most intense feature and accounts for the 3000A 
instrument efficiency cutoff, where appropriate. It can be seen in Table 6.01 that the 
spectrum is quite complex. Analysis of experimental results will require unfolding each 
constituent spectrum from the combined observed spectra. 
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TABLE 6.01 



MUSTANG EXPERIMENT; OBSERVABLE EMISSIONS 



Soecies 


System 


Band 


uAi 


Observable 


n 2 


VK 


(5,1) 


1841 


above 200km 




2P 


(2,0) 


2973 


marginal due low Q-o) 






(1,0) 


3157 


marginal due to inst. eff 




LBH 


(3,9) 


1854 


yes 


O 2 


SR 


(3,7) 


2442 


yes 




H-I 


(12,1) 


2528 


marginal 


NO 


Y 


(1,0) 


2149 


yes 




e 


(0,1) 


1950 


yes 




5 


(0,1) 


1980 


marginal 


O 


line 


N/A 


2972 


above 175km 



B. TOPICS FOR FURTHER INVESTIGATION 

This preliminary identification of the predicted experimental spectrum for the 
MUSTANG experiment must initially be confirmed by data analysis to identify those 
emission features in Table 6.01. This data will then be analyzed to identify N 2 density in 
the E, FI, and F2 regions along with NO density in the E region. Further analysis along 
with data from the HIRAAS instrument will be used to identify the solar EUV flux and 
photoelectron flux. These results will contribute significantly towards meeting the DOD 
goal of identifying the ionospheric electron density. The experimental wavelength range 
allows observation of all major neutral contributors and a self consistent model of the 
electron density may be possible when analysis of all experimental results is complete. 
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